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Abstract 


We have performed an exact numerical calculation of the linear growth 
and phase velocity of Kelvin-Helmholtz unstable wave modes on a supersonic 
jet of cylindrical cross section. An expression for the maximally unstable 
wavenumber of each wave mode is found. Provided a sharp velocity discontinuity 
exists all wave modes are unstable. A combination of rapid jet expansion 
and velocity shear across a jet can effectively stabilize all wave modes. 

The more likely case of slow jet expansion and of velocity shear at the jet 
surface allows wave modes with maximally unstable wavelength longer than or 
on the order of the jet radius to grow. The relative energy in different wave 
modes and effect on the jet is investigated. Energy input into a jet resulting 
from surface instability is discussed. 
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I. Introduction 


It has been suggested that surface instabilities such as the Kelvin- 
Helmholtz instability generate internal turbulence in jets and radio sources 
which can serve to reaccelerate particles (Pacholczylc and Scott 1976) and 
Dobrowolny (1972) has shown that the Kelvin-Helmholtz instability generates 
MHD waves in a collisionless plasma. Several authors have assumed the 
existence of turbulence in the form of' small amplitude MHD waves and used 
this to investigate particle acceleration. The wave spectrum was assumed 
to be a power law (Lacombe 1977; Ferrari, Trussoni and Zaninetti 1979) or 
related to the growth rate of surface instabilities (Eilek 1979). The MHD 
waves were found to heat the fluid and accelerate particles but did not 
lead to an appropriate particle energy distribution or emission spectrum 
if radiative losses were Included. These calculations also assumed a uniform 
input of wave energy throughout a source but energy input at the surface 
implies that wave energy is not uniformly distributed. This last problem 
has been addressed by Eilek (1982) who finds that energy input restricted 
to a thin surface layer produces a turbulent edge in MHD waves and strong 
limb brightening can occur. However, observed radio jets do not appear to 
be strongly limb brightened. A partial resolution of these problems may be 
found through investigation of the behavior of perturbations to a jet 
surface which is Kelvin-Helmholtz unstable. 

Supersonic jets in pressure balance with an external medium are unstable 
to surface perturbations analogous to the Kelvin-Helmholtz instability of 
a vortex sheet (Ferrari, Trussoni and Zaninetti 1978; Hardee 1979; Ray 1981). 
Growing perturbations of long maximally unstable wavelength are pinching 
and helical wave modes and previous wrk has provided estimates of wave phase 
velocity and maximum rate of growth (Ray 1981; Hardee 1982). In addition to 


XIX-13 


these wave modes a cylindrical jet with sharp velocity discontinuity is 
unstable to an infinite number of other wave modes with shorter maximally 
unstable wavelength and more rapid rate of growth. These growing waves 
convert directed flow energy into wave energy in the jet and external 
medium. Because initial growth of surface perturbations is rapid several 
authors have investigated ways of slowing the rate of growth and stabilizing 
a jet. The maximum rate of growth is decreased as the amplitude of a 
helical wave increases (Benford 1981) and growth of perturbations to the jet 
surface can be linear rather than exponential if a jet expands' rapidly 
enough (Hardee 1982).. Short wavelength perturbations to a jet surface can 
be stabilized by the presence of a velocity shear between jet and external 
medium (Ray 1982). While it may be possible to completely stabilize a jet 
by a combination of jet expansion and velocity shear, true jet conditions 
are not yet well known and it seems likely that some instability will be 
present. For this reason we must consider the behavior of long and short 
wavelength perturbations to a jet and consider the effect of the resulting 
growing wave modes on a jet. 

In this paper we will proceed by calculating exactly the wave phase 
velocity, growth rate and maximally unstable wavelength of Kelvin-Helmholtz 
unstable eigenmodes on a jet of cylindrical cross section. Both jet and 
external medium are assumed to be isothermal (see Hardee 1982, hereafter H II). 
The estimates made in H II forplnching and helical wave modes are confirmed 
and the work is extended to the wave modes describing perturbations of 
shorter wavelength. Stabilization of these wave modes is investigated. The 
depth interior to the jet surface affected by a wave mode is derived and 
relative energy input into a jet is estimated. The consequences for particle 
acceleration models are discussed. 
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I 


II. Dispersion Relation and Eigennodes 

It was shown in HII that the wave modes which propagate along the 
surface of a jet of cylindrical cross section are eigenmodes of the dis- 
persion relation 




In equation (1) is the Bessel function and is the Hankel function 

describing an outward propagating disturbance. The primes indicate derivatives 
with respect to the arguments of the Bessel and Hankel functions and the meaning 
of the symbols is 

E E y[(4> - - (1 - 4^ + 21(1 - 4^ O/yk]^^^ 


^n = - (1 + n^/k^ - 2i/k)v]^^^/v^^^ 

where 


4> = **^/ka . 


M * u /a . 
in in 


Y - (1-u^/c^)"^^^ 


'’ex * 


V - a ^ /a^. 
ex in 


The subscript in or ex refers to quantities interior or exterior to the jet 

respectively. This dispersion relation is valid for jets of constant radius 

or expanding jets of instantaneous radius R ■> r Y where r is the distance 

from the origin and ¥ = sin G where G is vone half the jet opening angle. 

In this notation k¥ ■ k R is the wavenumber normalized to the Instantaneous 

ex 

jet radius and k^^* 2x/A is the wavenumber. Eigenmodes of the dispersion 
relation are wave modes with n*0, 1, 2... . A wave mode propagates along 


the jet surface with form A ■ A exp (i(k r + n4> ” ut)]. 

n «x 
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The n * 0 wave mode pinches a jet and the n * 1 wave node helically distorts a 
jet. In HII it was shown that these two wave nodes were naxinally unstable at 
wavelengths <v, o.6(M^^/ri^^^)R and X^^^ 1.6(M^jj/n^^^)R. The wave node charac- 
terized by n ■ 2 distorts a jet in an elliptical shape; the wave node charac- 
terized by n * 3 distorts a jet in a triangular shape, etc.. For convenience we 
will call all wave nodes with n > 1, fluting wave nodes. These waves will 
not alter the jet cross section appreciably like the n * 0 wave node or nove the 
jet axis transversely like the n = 1 wave node. 

We expect fluting wave nodes to be naxinally unstable at sone wavelength 

X^^Nust as we found for the n = 0 and n = 1 wave nodes. To see that this 
n ■' 

nust be the case we take equation 1 and write it in the following fom 
(Appendix B of HII) 






^ -Vi (-n„)/-i„(kn„)i 

'-il; -"“i” 


( 2 ) 


When n/k'}'^ » 1 and n/kW^ » 1 equation (2) is easily seen to become 

n n 


_ ( . ) 

2 _ in 

Y n 2 — 


(3) 


This expression is just the usual expression for Kelvin-Helmholtz instability 
of a vortex sheet: 


0 ) ' V . [ 


y n 

i*hr^n 


± i 


Yn 


1/2 


1+Y^n 


]ku. 


(4) 
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We can see that the small wavenumber approximation used to obtain equation (4) 
remains valid to larger wavenumbers as n is increased. Thus we expect equation 
(4) to adequately describe the growth and propagation of waves with higher n 
to larger wavenumber, and we expect faster growth at shorter wavelength as 
n is increased. We can show that a maximally unstable wavenumber must exist 
by examining equation (1) in the limit n/k'J'^^«l .and n/kH'C^« 1. In this 
limit equation (1) can be written (Appendix B of HII) in the following form 





where Xj^= kV^^j ^ 

If we take the natural logarithm of equation (5) then 


V 1 


2n+l, 
4 ^ 


tr 

kV 


2k}¥ * 


(5) 


(6) 


where A^”^ is the logarithm of the right hand side of equation (5) . If equa- 


tion (6) is squared and we use and 


. . (n) . (n) 

assume that A^^ » A^ 


where and Aj"^are the real and imaginary points of A^”\ respectively, it 


follows that 




(2nfl)TT . (n) 

8k2y2 ^ 


( 7 ) 


We need to point out an error in equations (14) and (22) in HII in which the 
imaginary part needs to be multiplied by a factor of 2* For »R 1 ” 
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find that 0 and all wave nodes are unstable but with rapidly decreasing 

growth rate as the wavenumber becomes larger. Thus a maximally unstable 
wavenumber must exist. At a given wavenumber we see that the growth rate in 
equation (7) becomes larger as n increases. 

The maximally unstable wavenumber, maximum growth rate and wave phase 

velocity of the different wave nodes cannot be obtained analytically. We 

have found exact solutions to equation (1) numerically. In figure 1 we 

plot the growth rate in cm sec \ where R is instantaneous jet radius as 

a function of the normalized waven\imber kY « 2 itR/X for wave modes n equal 0 

through 5 and 7 and 9. The calculations assume a jet 10 times denser than 

8 “*1 

the surrounding medium, n ■ 10, that the jet velocity u»9xl0cms, 
that M. = u/a. ■ 9.9, and that the jet is expanding with an opening angle 
2G « 0.08 radian. We find that the maximally unstable wavenumber .is not 
sensitive to jet opening angle. These exact ' results when combined with the 
density and Mach number dependence of the maximally unstable wavenumber from 


HII which is valid when > 1 and M, » 1 show that 

in 

k^P^ y '^1.9(n+ -377^^= 

n + 0.2 


At large n^k^^^'F 'v- 1.9n(Ti^'^^/M^^) . The asymptotic expression is accurate to 
10% when n=3. 

The corresponding relationship for the wavelength is 


3.3 [ 


■■■ 3/2 1® 

n+l.l/(n-*'^-K).2) 


For n less than 3 we have that 0.6 (M£^/n^^^)R, ^ 

and fv 1 .4 (Mj^/ti^^^)R. For larger values of n,X^^^'^3.3n is 

accurate to better than 10%. The values of X^^^ and x5^^ found here are in 

O X 

excellent agreement with the estimates given in HII. 
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The phase velocity of the different wave inodes noznallzed to the jet 

velocity u as a function of the normalized wavenumber is shown in figure 2. 

At wavenumbers much less than the maximally unstable wavenumber the phase 

2 2 

velocity, v^j^, approaches u or [y r)/(l+Y .ti]u for the n • 0 and n > 0 wave modes, 

respectively. At wavenumbers much larger than the maximally unstable wavenumber 

the phase velocity approaches u - a. for all n » 0. For the n = 0 wave mode 

in * 

the minimum phase velocity occurs at wavenumber larger than the maximally 

unstable wavenumber. The phase velocity at the maximally unstable wavenumber 

is Vpj^ “v u-a^^. This is the result for a small amplitude wave driven by the 

Bernoulli effect. Larger amplitudes will develop only for smaller wavenumbers 

and higher phase velocities and any large scale jet pinching should develop 

with wavelength greater than • For all other wave modes the' minimum 

phase velocity occurs near the maximally unstable wavenumber. In general, 

when n > 1 we may assume that the dominant wavenumber and wavelength associated 

with a particular wave mode is within 10% of the value given by equations 

(8) and (9). We note that the reduction in phase velocity seen at the dominant 

wavenumber implies more rapid spatial growth than if the phase velocity re- 
2 2 

mained between [y q/ (l+y n] u and u-a. as assumed in HII. For n ■ 1 the 

in 

minimum phase velocity is about 10% less than the value predicted using equation 
(4). This suggests that spatial growth of the helical wave is 10% faster than 
was estimated in HII. For other wave modes as n becomes larger the maximum 
growth rate approaches that predicted using equation (4) with k ■ 

Vp^ is reduced significantly below that predicted using equation (4) with result that 
spatial growth is somewhat more rapid than would be predicted using equation (4). 


III. 


While this result depends on M. and h we find no large change as M. 

in in 

and ri are varied. As a consequence when n is sufficiently large, we expect 
initial wave node growth at small wave amplitude to be approximated by equation 
(4) until non-linear effects at larger amplitude become important. Of course, 
if jet expansion is rapid enough the wave remains in the small amplitude limit 
and spatial growth is linear as the maximally unstable wavelength increases 
proportional to R = rV (see HII) . 


Jet Application 

Jet expansion' can serve to stabilize the wave nodes which are maximally 
unstable at longer wavelength- To see that this cannot be the case for wave 
modes that are maximally unstable at shorter wavelength recall that 
[Y^o/( 1+ Y^n)]u and Uj < [Yh^^/(1+Y^h) ] 2 jtu/X^^^ for all wave modes with n > 1 . 
We define a spatial growth length as * v where 


equation (4) exactly to obtain v^^ and and -use 3.3 n ”')R 

then £ ^<0.52 n ^ M. R. This value of the spatial growth length is suffi- 

e in 

clently accurate for our purposes and exact calculation shows the numerical factor 
is only a weak function of n. Wave amplitude grows as A ® A^e "where r 

is distance along the jet surface. For an expanding jet R ® rV and (r/£^) “v. 

1.9 n/(r)^'^^M. 'f) . While the growth rate and phase velocity of growing pertur- 
batlons will change as wave amplitude increases it is certain that a large 
value of (r/£^) implies that large wave amplitudes can occur. For our par- 
ticular choice of parameters 'v- 0.85 and (r/£^) 2.2 n. It is 

clear that the helical wave mode (n ■ 1) can be restricted to only a few 
e-folding lengths but this is not the case for n > 1. While somewhat more 
rapid jet expansion is of some help we point out that < 1 for a jet in 


- 1 , 


-l/3> 
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I 

,1 





pressure balance with the external medium. A Jet in free expansion with 

1 is stabilized by virtue of its expansion at the sound speed which 

restricts surface instability to a thin surface region. 

If a jet is in pressure balance with an external medium the only vay 

to stabilize the jet surface to wave modes with shorter wavelengths is by 

velocity and density shear . For a linear velocity shear layer Ray (1982) 

finds that waves with wavenumber k > 1.28 h ^ where h is the scale height 

* 

of the shear layer are stabilized. To apply this result to a cylindrical 

jet we need to determine the true wavelength of a wave mode at maximiim growth. 

The true wavelength is always less than the wavelength X at constant angle 

(J) (cylindrical or spherical geometry) . Figure 3 illustrates the propagation 

geometry of wave modes on a Cylindrical jet of constant cross section. The 

2 2 1/2 

true wavelength between wave fronts is given by /n]/l (2irR/n) +X^) ' )X^. 

The true maximally growing wavelength for a wave mode and therefore the true 
characteristic wavelength of a wave mode is 

2 1/2 

- {[2TTR/n)/I(2TTR/n)^ + ] }X^P^ . (10) 

n n . n 


This relation is modified only slightly by jet expansion (see figure 1 in HII) 
and may be used with confidence with R being the instantaneous jet radius. 

The stability criterion can be written A < 4.91 h. If* for example, the depth 
of a shear layer were R/4 then A < 1.2 R would be stabilized. It is informative 
to rewrite the stability condition in terms of X and we find that %»avelengths 

X > X“^" »2TTn”^ I(1.28R/nh)^ (11) 

n n 


. are unstable and shorter wavelengths are stabilized. Immediately we see that 
wave modes with n > 1.28 R/h must be completely stabilized. This occurs 


because no matter the wavelength X^ for these wave modes, 

A is less than the minimum unstable wavelength . 
n 


the true wavelength 
Since the maximum 
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depth of the shear layer is h *=^ R for a cylindrical jet it would not appear 
possible to stabilize the n « 1 helical wave mode by velocity shear. However, 
it is likely that only wave modes with small n, say n < 10, will be unstable. 
For example, using the parameters f) » 10, = 9.9 and h * R/4 we find 

7.8R > X®^” - 3.1R, through X^**^ 3.8R > X®^” - 3. OR but X^^^-v 

3.0R < X®^" » 8.1R. The maximum rate of growth of the wave mode n » 5 is 
reduced to about that of the n = 2 wave ipode (see figure 1) and will actually 
be reduced to less than this if full effects of a shear layer are included 
(see Ray 1982). Wave modes with n > 5 would be absolutely stable. We conclude 
that only a few wave modes will be unstable on an expanding jet with signi— 


ficant velocity shear. 

The observational consequences of the n * 1 helical wave mode have 
already been discussed in detail in HXI. It is now necessary to consider 
the effect of the remaining unstable wave modes on the jet material* Tne 
effect of shorter wavelength wave modes on the jet material will be 
limited by the portion of a jet affected by a surface wave. For a jet of 
cylindrical cross section the displacement amplitude of jet material interior 


to the jet surface is given by (HII equation A18) 

2 


1 r ^in ^in 


n(e) - 7 I — 2 i 

“in 'I’in * > 


J IB*-” cose 


j' (B^” « 

n n 

j (b!;" 4') 

n n 


In equation 12 6^” H and spherical geometry is assumed. At the jet 
^ n n 

surface 6*0 and 4; * To obtain the amplitude as a function of jet radius 


we form the ratio 


n(6)/ri{0)> 


^ J • (bJ," 4^) 

cos 6 n n 

cos 0 j’ in^-) 

n n 


where Pi„ (6)/Pin(e) = 

obtain an estimate of the depth affected by each wave mode n. 


( 12 ) 


( 13 ) 
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Provided n > 2 the characteristic wavenuaber for each vave.node 
can be approximated by 'v l.9n with acceptable accuracy. 

The arguments of the Bessel functions in equation 13 are proportional to n 
and we may use the following approximation for the Bessel function: 

J„(nz) 'V' ^ Ai 

1-2 


( 14 ) 


where Ai is the Airy function and 


2 3/2 
3^ 


£n 


l+'^l-z^ 


- /1-z" 


or 


f /z^-1 cos (1/z) 


so X is real when z is positive (Abramowitz and Stegun 1964) . In equation 


,1/3, 


(14) 2 • 1.9 (ti /M. )(tf)/'F) If we use equation (4) for 4* * t- at k 

in n ’ ka. n 

in 

and the fact that 1 for a supersonic jet it follows that 


(P) 


L-fM 




1/2 


n « 1 
n » 1 


As a result we find 


2 ' V , 


rx., 

1,9 




^ n 


- 1/6 


n « 1 


/.I. /ti/\ 




and in general |z|<l. When |z| is small X ( - |ln(2 /z)-l) AKn^^X)*^ 

expl“-| nx^^^] and 4x/(l-2^)*'' 4x(l+z^) can be used with 

result that 


J (nz) *v — 
n r 


(«_) 


n 


ZTin 


( 15 ) 
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Because J ’ (nz) * Ij (nz)-J (nz) (nz)/z it follows that 
Ti n n-rl “ 

z 

t' r ^ (e/2)^ n-1 

The approximate fonn of the derivative of the Bessel function inserted 
in equation (13) results in the following relation between the displacement 
amplitude at the jet' surface and inside the jet: 

n-1 


cos 0 

h(e)/Tl (G) eos 


( 16 ) 


(17) 


Equation (17) implies that the wave amplitude interior to the jet surface is 


A -v (b/R)"~^ 

n n 


(18) 


where is the surface amplitude and 0 <b <R where b is the radial distance 

inside the jet. The lack, of radial dependence of the n * 1 helical wave mode 

results because the entire jet is displaced. 

If conditions are such that (r/£ ) >> 1 for a wave mode then it is 

e 

reasonable to assume that the wave amplitude at the jet surface has grown 

to its maximum. While maximum amplitudes and associated fluid motions need 

to be determined by numerical simulation, ordered displacement of the jet 

surface is limited to motion of the fluid at the sound speed transverse to 

the flow velocity. Wave modes with n > 1 appear as sinusoidal displacement 

of the jet surface r](r) = sin(27Tr/X^^^) at fixed angle • The maximum 

^ ^ R 1 ( ) ^i 

amplitude of such a sinusoidal oscillation is A :: — ] where 

n xtt n u — V , 


ph 


(p) 


u - v^^ is the fluid velocity in the wave frame. Because v^^ ^ n 

R -1 

only a weak function of n, A^ « n . The energy in each wave mode is 

2 

proportional to AM(Au) where AM is the jet material participating in motion 
at velocity Au. If we continue to assume that Au a^^ at the jet surface 


.n-1 


then Au « (b/R)** independent of To the extent that wave modes are 
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independent and all of the jet material is involved » the energy in a wave 

mode interior to the jet surface is E “ (b/R)^^*' The total wave 

n 

energy interior to the jet surface is of the form 

E 'v E y a (b/R)^^" " (19) 

o ^ n 
n 

A fundamental scale size of K exists for each term in the summation and 

n 

a is a weak function of k . An uppef limit to the summation is set by 
n n 

the depth of the shear layer and a lower limit may be set by jet expansion. 
Near the upper and lower limits must decrease rapidly. Near the upper 
limit wave growth is slowed because the fundamental scale size is shifted 
to longer wavelength and at the lower limit wave growth is slowed and 
maximum wave amplitude is not reached. 
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IV. Discussion 


Each unstable wave mode can be regarded as driving the jet surface at 
a length scale and a frequency • (u - 

is less than the jet diameter it is likely that eddies of size can 

be created in the jet fluid. Dissipation in the jet fluid leads to a 

“5/3 

cascade to eddies of smaller size with a Kolmogorov energy spectrum « k 

and with the smallest eddy size governed by dissipation processes in the 

fluid (Landau and Lifschitz 1959). Wave modes with larger than the 

jet diameter cannot generate eddies in the jet and must affect the flow in 

a more global sense. If we consider n * 10 * 9.9 we find that 

» 4.9 R, A 2 ^^ » 2.8 R and A^^j * 5.8 n“^R. In this case the helical 

and elliptical wave modes cannot produce eddies in the jet material without 

disrupting the flow- These wave modes must dissipate jet energy in the 

external medium and jet by heating of the fluid (Benford 1981). However, 

each unstable wave mode with n > 3 can product eddies in the jet material. 

Other values of ri and can reduce less than a jet diameter. 

To the extent that each wave mode is independent, each drives a fundamental 

eddy size with total energy associated with each wave mode given by the 

appropriate term in equation (19). Wave modes which produce eddies in the 

jet dissipate energy at a rate -v- (Au)^/A^^^ (Landau and Lifschitz 1959). The 

dissipation rate increases as the fundamental scale size decreases. 

If we continue to assume that Au ^ the jet surface independent of 

then energy input into the jet is of the form 
n 

J -V e J B (b/E)’‘" • ” • (20) 

o n 
n 

where 6 ^ n a . 

n n 
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In equation (20) a lower limit to n is set when is greater than the let 

tl 

diameter as dissipation is no longer through fluid turbulence driven at scale 
• (p) 

size A ^ . Modifications to this simple relation will occur if wave modes 

strongly interact and if fluid turbulence generated at shorter A^^^ propagates 

into the jet faster than perturbations at longer A^^^ grow to large amplitude. 

The effect of Kelvin-Helmholtz unstable surface waves on the jet interior 

and external medium is of fundamental importance to an understanding of 

observed jets. We have found that the presence of velocity shear may stabilize 

all but one or two of the long wavelength wave modes and that these long 

wavelength wave modes can be restlcted to only a few e-folding lengths and 

small amplitude by jet expansion. If these conditions are met then we might 

expect that little of the flow energy would be dissipated in the external 

medium or internally in the jet because turbulence would be unlikely to 

develop. However, if the jet edge is relatively sharp then we expect a 

number of wave modes with A^^^ less than the jet diameter and (r/£ ) >> 1. 

n e 

The likely case would be one in which some wave modes with say 2 < n $ 10 

grow rapidly and those of larger n are surpressed by velocity shear. 

It is clear that restrictions on wave mode stability are likely to result 

in fundamental eddy sizes on the order of the jet radius or suppress the 

development of short wavelength wave modes and jet turbulence entirely. If 

wave amplitude is limited by motion of the jet surface at the sound speed 

then wave energy is proportional to (Au)^ u^/M? and can be a large or 

in 

small fraction of the jet flow energy. Interestingly, if » 1, velocity 
fluctuation in the jet material at the sound speed is a small fraction of 
the flow speed and the flow can appear well ordered in the observer frame. 

Even so, the flow would be characterized as strongly turbulent if lOZ of 
the flow energy were converted to eddies in the jet frame. 
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Our present results suggest that surface Instability can affect more 

than a thin surface layer and probably produces a Kolmogorov energy spectrum 

in fluid turbulence in a 5^eti' Some of the wave energy converted into 

turbulence will simply heat the fluid through dissipation but the presence 

of magnetic field implies the generation of MHD waves (Kato 1968; 

Dobrowolny 1972). If the fluid is driven turbulent at a single wavenumber 

or only a narrow range of wavenumbers the MHD waves produced will have an 

-3/2 

energy spectrum “ h (Kraichnan 1965; DeYoung 1980). A Kolmogorov 
energy spectrum in fluid turbulence allows direct drive to MHD waves at 
many wavenumbers. The resulting MHD wave energy spectrum is not well 
determined but is probably as steep or steeper than Kolmogorov (Kato 1968; 
Henriksen and Eilek 1982). Whatever the exact spectrum, the MHD waves 
produced in this fashion no longer are confined to a thin surface region 
(see Eilek 1982) but can exist in a significant fraction of the jet volume 
with MHD wave energy proportional to the wave energy in fluid turbulence 
(see equation 19) . This MHD turbulence can supply energy to an existing 
population of relativistic electrons and may produce the observed jet 
emission (Eilek and Henriksen 1981) . Since MHD waves of the appropriate 
wavelengths for particle acceleration are strongly damped and do not 
travel far (Eilek 1982) the jet emissivity may be of the same form as 
equation (20). Such a jet emissivity should result in jets with intensity 
profiles flatter than if the emissivity were not a function of the jet radius. 
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2i 

Figure 1: Wave mode growth rate (jOjR in cm sec as a function of vavenuisber 

normalized to the jet radius on a jet 10 times denser than the surrounding 
medium, Mach number of 9.9 and opening angle 20 « 0.08 radian. 
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Figure 3. Propagation geometry of wave fronts (dashed lines) of wavemode^ 
n = 2 and 3 around a cylindrical jet of circumference 27 tR at equal propaga- 
tion angle with respect to the flow velocity. 


XIX-33 


